- 1 - 



Dynamical magnetic relaxation: A nonlinear magnetically 

driven dynamo 



Eric G. Blackman 1 and George B. Field 2 
1. Department of Physics & Astronomy and Laboratory for Laser Energetics, University of 

Rochester, Rochester NY 14627 
2. Harvard-Smithsonian Center for Astrophysics, 60 Garden St. Cambridge MA, 02138 

(in press, Physics of Plasmas) 
ABSTRACT 

A non-linear, time-dependent, magnetically driven dynamo theory which shows how 
magnetically dominated configurations can relax to become helical on the largest scale 
available is presented. Coupled time-dependent differential equations for large scale 
magnetic helicity, small scale magnetic helicity, velocity, and the electromotive force 
are solved. The magnetic helicity on small scales relaxes to drive significant large 
scale helical field growth on dynamical (Alfven crossing) time scales, independent of 
the magnitude of finite microphysical transport coefficients, after which the growing 
kinetic helicity slows the growth to a viscously limited pace. This magnetically driven 
dynamo complements the nonlinear kinetic helicity driven dynamo; for the latter, the 
growing magnetic helicity fluctuations suppress, rather than drive, large scale magnetic 
helicity growth. A unified set of equations accommodates both types of dynamos. 

PACS codes: 52.30.Cv, 95.30.Qd, 52.65Kj, 52.55 -s, 96.60.Hv, 98.62Mw 



I. Introduction 



The concept of dynamo means different things to different communities: In the astrophysical 
context, dynamo typically refers to the amplification of initially weak magnetic fields in a turbulent 
flow, but comes in two basic varieties: (1) small scale nonhelical (direct) dynamos in which 
magnetic energy is amplified by random walk field line stretching on spatial or temporal scales less 
than or equal to that of the turbulent forcing [1-4], and (2) large scale helical (inverse) dynamos 
[5-12], in which helical turbulence produces magnetic energy on scales larger than that of the 
input turbulent forcing. The latter is the type of dynamo needed to explain the large scale field 
and solar cycle of the sun. It thrives on a non-dissipative term in Ohm's law called the turbulent 
electromotive force (EMF), E , which equals the correlated average (time, space or ensemble) of 
the cross product of fluctuating velocity v and fluctuating magnetic field b, and is sustained by 
helical turbulence. 
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In fusion devices, the plasma is already magnetically dominated. Here the dynamo acts to 
convert magnetic flux from toroidal to poloidal (or vice versa), increase the scale of the field, and 
sustain the strong magnetic flux against microphysical dissipation [13-21]. Like the velocity-driven 
astrophysical helical dynamo, the magnetically dominated dynamo thrives on having a finite £ 
term in Ohm's law, again resulting from correlated turbulent fluctuations possessing helicity (more 
on this later). But here the fluctuations are driven by current instabilities and magnetic helicity 
injection rather than kinetic helicity. 

Although traditional astrophysical dynamos are normally studied for systems with initially 
weak magnetic fields, there are also important magnetically dominated astrophysical environments: 
astrophysical coronae. These include solar and stellar coronae, as well as coronae above accretion 
disks [23-25]. In coronae, the dynamo can convert and sustain flux, playing a similar role to the 
magnetically dominated dynamo in fusion devices. Coronal systems are magnetically driven in the 
sense that helical magnetic fields are injected from the disk or star below, and the subsequent time 
evolution is of interest. A key question is: how does the field open up and relax to large scales? 
This is relevant for the formation of magnetically driven jets and outflows from disks and stars. 

It is useful to recognize that a magnetically dominated dynamo can also be thought of as 
dynamical magnetic relaxation. Magnetic relaxation describes the process by which magnetic 
structures in magnetically dominated environments evolve to their equilibrium states. The fully 
relaxed end state is the Taylor state [22], determined by minimizing the magnetic energy subject 
to the constraint that magnetic helicity is conserved. The result is a force-free helical configuration 
with the scale of the field reaching the largest scale available, subject to boundary conditions. 
But Taylor's theory by itself is not a dynamical theory since it does not provide a time-dependent 
description of how the large scale magnetic helicity evolves. For that, a fully time dependent 
dynamo theory is needed. That is our goal herein. 

In Ref. [12], a set of nonlinear dynamo equations for the growth of the large scale magnetic 
field was derived and solved when helical turbulent velocity forcing is applied. The nonlinear 
backreaction due to the build up of small scale magnetic helicity was shown to ultimately quench 
the large scale field growth, in quantitative agreement with numerical simulations [10, 11]. In the 
present paper, we study the complementary problem of large scale field growth from a dynamo 
driven by small scale magnetic helicity fluctuations. As we will show, the complementarity arises 
because the kinetic helicity becomes the quenching agent rather than the driving agent for the 
magnetic helicity driven dynamo. A novel feature of the present approach is that we include a 
fully time-dependent dynamical equation for the turbulent EMF which couples magnetic helicity 
and kinetic helicity, dynamically into the time dependent theory. Previous work on magnetically 
dominated dynamos focused on more specific magnetized configurations [13-21] in the steady-state 
in which S is not solved for dynamically. We find that coupling the time evolution of £ into the 
theory is essential for understanding time-dependent magnetic relaxation. 

We solve the dynamical nonlinear dynamo equations for a closed or periodic system. This 
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allows us to ignore boundary terms and will facilitate future testing of the basic principles with 
tractable 3-D numerical simulations. The theory is subtle enough, with enough new features, that 
we wish to present it for as simple a system as possible without focusing on detailed magnetic 
configurations and boundary terms. These will have to be considered in future work. 

In section II we derive the coupled equations to be solved: the time evolution of large and 
small scale magnetic helicities, the turbulent EMF, and the kinetic helicity. In section III we 
discuss the solutions and the physical interpretation, and conclude in section IV. 



II. Derivation of the Coupled Dynamical Equations to be Solved 

We will use a simple multi-scale approach that has been relatively successful [12] in accounting 
for the non- linear dynamics of recent 3-D MHD helical dynamo simulations [10, 11] (and has 
been further been shown to be reasonably consistent even compared to a theory in which 2 
additional scales are included in the dynamics [26]): We write all quantities as the sum of large 
scale (indicated by overbar) and fluctuating (indicated by lower case) contributions. The vector 
potential A, magnetic field B = VxA (written in Alfven units) and normalized current density 
J = VxB then satisfy A = A + a, B = B + b and J = J + j respectively, where the overbar 
represents a local spatial average. We assume a separation of scales such that the lower case 
quantities vary on scale of inverse wavenumber k^ 1 and the overbarred quantities B, A vary on 
scale k^ 1 (0 < k\ <C fo). We also define a global spatial average taken over the system scale 
(3> k^ 1 ), or periodic box. The global averages are indicated by brackets (), and being global 
averages, have no spatial dependence in our approach, only a time dependence. 

We define the total average magnetic helicity H M = (A • B), and large and small scale 
averaged magnetic helicities as = (A ■ B) and = (a • b), such that H M = Hf + Rf. 
The induction equation 

8 t B = -VxE = Vx(v x B) + AV 2 B, (1) 

(where v is the fluctuating velocity, A is the magnetic diffusivity) and E = — <9 t A — V^>, imply 
that [6] 

d t H M = dtHf* + dtHf = -2(E • B) = -2A(J • B), (2) 

where <fi is the scalar potential. The large and small scale integrated magnetic helicity equations 
satisfy 

d t H^ = 2(£ ■ B) — 2A(B • J), (3) 

and 

d t H^ = -2(S -B)-2A(b.j), (4) 

where the turbulent electromotive force £ = v x b = — E + A(J • B). Note that in deriving (4), we 
have used (2), (3), and the definition of £. 
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For later use, it is helpful to derive relations between globally averaged current helicity, 
globally averaged magnetic helicity, and globally averaged magnetic energy. For the current 
helicity we have 

(J • B) = -(B • V 2 A) + (B • V(V ■ A)) = -(B • V 2 A) (5) 

where the last equality follow from the chain rule, V • B = and the fact that total 
divergences vanish when globally averaged. Correspondingly, |(b • j)| = \k\H^f \ < ^(b 2 ) and 
|(B • J) | = |/c 2 i?f f | < fci(B ). The inequalities in the preceding two relations represent the the 
realizability condition [27] and the equality holds for maximally helical (force-free) structures on 
the respective scales. (These are gauge independent relations because globally averaged divergences 
vanish. If we had used the Coulomb gauge (V • A = 0), then (5) and the relations above would 
have followed straight away, without having to appeal to vanishing global divergences.) 

The dynamical equation for the turbulent £ is 

d t £ = dtv x b + v x 8 t b. (6) 

Therefore, we need equations for dth and div. Assuming that V • v = 0, subtracting the local 
mean of (1) from (1) we have 

d t h = Vx(v x B) + Vx(v x b) - Vx(voTb) + AV 2 b. (7) 

Similarly, in the absence of any mean velocity, the equation for fluctuating velocity becomes 

d t v = v x w-vxw - V(p + v 2 /2 - v2/2) + j xB + Jx b + iA7 2 v, (8) 

where v is the viscosity and p is the fluctuating pressure. In deriving (8) we assume that j x b = 0, 
from being injected with this property. For simplicity, we also assume that this is maintained at 
all times (this needs to be tested). This assumption should not affect the basic conclusions from 
our study because only j • b enters the turbulent electromotive force £. In addition, we do not 
assume that the magnetic force associated with the cross-terms between large and small scales are 
force-free, and thus we still allow terms such as j x B and J x b to contribute to the small scale 
magnetic force. 

By analogy to Ref. [12], using (7) and (8) we then have for (6) 

d t £ = i(b • Vxb - v ■ Vxv)B - pVxB + vV 2 \ x b + Av x V 2 b + T, (9) 

where Tj = [v x Vx(v x b)]j — ej qn P q i{v ■ Vvi)b n represents the surviving triple correlations (by 
assuming j x b = there no triple correlations survive with more than one factor of b) and where 
P q i = (S q i — V~ 2 V g Vj) is the projection operator that arises after taking the divergence of the 
incompressible Navier-Stokes equation to eliminate the pressure. In deriving (9) we took v and b 
to be statistically isotropic. It will be of interest in future work to drop this assumption [28, 29]. 

For Eqns. (3) and (4), we want the component of £ parallel to B. For this we have 

8 t £\\ = (ftvxb + vx d t b) • B/|B| + Wb • d t (B/\B\). (10) 
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Using Eqns. (9) and (10) we have 

d t £\\ = dB 2 /|B| - (3B • VxB/|B| - (£ „ (11) 

where a = (l/3)(b • Vxb — v • Vxv), f3 = (l/3)v 2 , and C = /^v 2 *^ accounts for microphysical 
dissipation terms, the last term of (10), and most importantly, T. We will take the constant 
/ ~ 1, which follows from estimating the magnitude of T. Note that £ and (3 depend on v 2 which 
must be solved for dynamically. 

When B is force-free and only one sign of helicity is initially injected into the volume, 
B • J ~ (B • J) and B 2 ~ (B 2 ) = ki\H^\ [10]. (We consider the injection of one sign to illustrate 
as simply as possible the tendency for magnetic helicity to inverse cascade. If both signs are 
injected with equal magnitude then no net magnetic helicity would be injected at all. Also, the 
source of large scale field growth in the problem we consider comes from the £ terms, which grows 
force-free large scale fields.) In addition, v 2 ~ (v 2 ), v • Vxv ~ (v • Vxv), and b • j ~ (b • j). Then 
(3) and (4) become 

d t Rf = 2£\\k\ /2 \H^\ l l 2 - 2\k 2 H^ (12) 

and 

d t H™ = -2f||fc 1 1/2 |i7 1 M | 1 / 2 - 2\k\Bf. (13) 
Using U = (v 2 ), Hrf = (v • Vxv), we re- write (11) in the two scale approach as 

dtS\\ = k\ /2 \Hf\^ 2 (k 2 H^ - F 2 y )/3 - fc? /2 (fff /(|fffY/2)i7/3 - fhU 1 / 2 ^. (14) 

To solve for the evolution of and H^ 1 using (14), we also need dynamical equations for U and 
H\ . Our substitution v 2 ~ (v 2 ) allows us to ignore boundary terms, since we have stated that () 
is taken over a closed or periodic volume. Then, from (8), 

t U = 2((v x j) ■ B) — 2<(v x b) • J) - 2^((Vv) 2 ) ~ 2{k 2 - ki)(£ ■ B) - 2vk\(v 2 ), (15) 

and 

d t H\ = 2B • {Co x b)(fc 2 - h) ~ 2v{d i v j d i uo j ) ~ 2k 2 {k 2 - fci)(5 • B) - 2u{d i v j d i uj j ), (16) 

where we consider the case that Hf 1 , > 0, took advantage ofjxb = JxB = 0, and again 
dropped surface terms. We also used dt(v ■ h) = —k 2 {\ + ■ b) (showing that (v • b) decays) to 
write (Co x h) q = (v s d q b s ) = (v x j) 9 in deriving the penultimate term of (16). 

The equations to be solved are thus (12), (13), (14), (15) and (16) after converting them 
into dimensionless form. We write Hrf at t = as 11^(0) and define the dimensionless 
quantities h x = H^/H^(0), h 2 = i/ 2 M /i7 2 M (0), R M = ^Hplfi) / \k] /2 , R v = ^// 2 M (0)/^ 2 1/2 , 

r = tkl /2 ^/H 2 M {0), Q = -£ l] /k 2 H 2 M (0), e = U/hH^iO), and h v = i7 2 y /fc^ 2 M (0). Using these in 
(12), (13), (14), (15) and (16) respectively gives 

3M = -2QhT(kilk 2 f t2 - 2h 1 (k 1 /k 2 ) 2 /R M , (17) 
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8 T h 2 = 2Qh\ /2 {k l /k 2 ) 1/2 - 2h 2 /R M , 
d T Q = ~ (h/k 2 ) 1/2 h\ /2 (l/3)(h 2 - h v ) + (h/k 2 )V 2 h{ /2 e/3 - e^fQ, 
8 T e = -2(1 - ki/k 2 )(ki/k 2 ) 1/2 Qh\ /2 - 2e/R v , 



(18) 
(19) 
(20) 



and 



d T h, 



2(1 



1/2 



2h v /R v 



(21) 



i 



(In the simple treatment above we have not coupled in an equation for the non-helical part of the 
magnetic energy. This would arise for example, from a bulk large scale shear, which we also do 
not include. A varying fraction of non-helical magnetic energy on the large scale would add a time 
dependent coefficient > 1 to Q in (17) [30].) 

Below we consider two cases CI and C2: For CI, h 2 takes on an initial value h 2 (0) = 1, and 
then evolves according to (18) while for C2, h 2 = 1 is fixed for all times. Case CI corresponds to 
free relaxation (simulated in Ref. [31]) and case C2 corresponds to driven relaxation. The latter 
represents injecting h 2 into the system, for example by a potential difference [20, 21] which drives 
instabilities that keep h 2 steady. In both cases e and h v satisfy the same equation, so h v = e for 
all times when both are initially small. 



Solutions for case CI with / = 1 and different Rm = Rv are shown in Fig. 1. We see that hi 
grows and h 2 depletes in accordance with magnetic helicity conservation, while at the same time 
e = h v grows. This prevents h\ from equaling h 2 (0) = 1 due to the backreaction from h v in (19). 
Instead, the maximum hi occurs at h\ ~ h 2 ~ 1/2, after which resistive terms take over and all 
quantities decay. If h v did not grow, then h\ would reach h\ ~ h 2 (0) = 1. 

Solutions for case C2 with / = 1 are shown in Fig 2. In Fig. 2a and 2b, a significant 
kinematic regime of duration r^in in which growth of h\ is independent of Rm and Ry is evident. 
The quantity Ttin can be estimated by ignoring the Rm and Ry terms and replacing the time 
derivatives in (17), (19), (20), and (21) by multiplication by 1/rkin- The end of the kinematic 
phase occurs when h v grows large enough to deplete the second term in (19) such that the third 
term in (19) becomes important. This occurs when h v = e = (1 + ki/k 2 ) _1 . Using this, and 
solving for r^n from the approximated versions of (17), (19) and (20) gives 



which is ~ 40 for k 2 jk\ = 5 and / = 1. Solving for h\ at t^m, then gives hi ~ (1 — k 2 /k 2 ) 1 ~ 1. 
Thus hi grows kinematically up to a value ~ h 2 = 1. 

After Tfcj n , Rm and Ry become important. The ratio determines the saturation value 
of hi, while Ry determines the rate of approach to the saturation. That Ry rather than 



III. Discussion of Solutions and Physical Interpretation 



r fcm ^(3//2)(fc 2 Ai) 2 (l + ^iA2) 1/2 



(22) 
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Rm determines this rate results because the build up of h v suppresses growth of hi and Ry 
appears in (21) in the diffusion term. We can estimate the final saturation values expected 
for case C2 by setting the right sides of (17), (19), (20), and (21) equal to zero. The 
solution is e = h v = (1 + k 1 /k 2 )~ 1 , Q ~ -(k 1 /k 2 ) 1/2 (l - k 2 ilkly l l 2 (RyR M Y xl2 < 1, and 
h\ = (Rm / Rv){k2 1 k\) 2 {\ — kf/k 2 )^ 1 ~ {Rm / Rv){k 2 /ki) 2 . The saturation value of hi can be seen 
to depend on Ry /Rm in Fig 2c. 

For both cases CI and C2 we have also solved the dynamical equations for / = 1/10 (C2 
solution plotted in Fig. 2d). This corresponds to a closure in which the damping time in the last 
term of (19) is ~ 10. Oscillations appear at early times, similar to what is seen in Ref. [12] for the 
h v driven dynamo when / < 1 and when the closure is applied to passive transport of a scalar 
[32, 33]. Unlike the case of Ref. [12] however, h v evolves here, so the similar sensitivity to / is 
noteworthy. Here, for / < 1, Q oscillates around zero because h v can grow to be larger than hi. 
The sign of d T Q then changes and eventually Q changes sign as well. For / < 1, the last term 
in (19) cannot damp growth of Q fast enough to prevent its sign change. Eventually, the Ry 
term takes over in (20) and (21), damping the growth of e and the oscillations. The oscillation 
amplitudes provide a prediction for a given / and a diagnostic for the closure. Note that the rise 
to the first peak in hi corresponds to the end of the kinematic regime. Eqn. (22) shows that 
Tkin decreases with /, explaining why the first peak occurs earlier for / < 1 than the end of the 
kinematic regime for / = 1 (compare Fig. 2d with Fig. 2b). 

An important feature of the system of equations (16-20) is that they represent a unified 
framework for understanding aspects of BOTH kinetic helicity driven helical dynamos and 
dynamical magnetic relaxation depending on whether the source driving the fluctuations initially 
supplies kinetic helicity or magnetic helicity. If kinetic helicity is steadily supplied, we have the 
nonlinear dynamo of Ref. [12]. Then dfh 2 = dte = and Q is initially dominated by the h v 
contribution in (18) which drives growth of the large scale magnetic helicity in (19). The near 
conservation of magnetic helicity then means that the small scale magnetic helicity h 2 must grow 
with opposite sign to compensate the growth of hi. The growth of h 2 , in turn, offsets the kinetic 
helicity in the first term on the right of (18), slowing the growth of hi in (17), and quenching 
the system into a steady-state. In contrast, as shown above, when the driving initially supplies 
h 2 instead of h v , Q is initially dominated by h 2 . Then h v grows, offsets h 2 in (18) which in turn 
quenches the growth of hi and again drives the system into a steady-state. In this sense, the kinetic 
helicity driven dynamo and the magnetically helicity driven dynamos are fully complementary, 
and their non-linear evolution is explicable by the same set of equations with different initial 
conditions. In both the kinetic and magnetic helicity driven dynamos, the large scale magnetic 
helicity incurs an initially fast growth phase (kinematic regime) which is independent of the value 
of the finite microphysical diffusivities, as long as the driving fluctuation scale is far enough above 
the diffusive scales. 



One important difference between two types of dynamos is that for the kinetic helicity 
driven dynamo, the large and small scale magnetic helicities have opposite signs, whereas for the 
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magnetically driven dynamo, the large and small scale magnetic helicities have the same sign. 
This is easy to understand and highlights the role of magnetic helicity conservation: In the kinetic 
helicity driven dynamo, there is initially negligible total magnetic helicity, and the kinetic helicity 
acts to drive one sign to large scales and the other sign to small scales. In the magnetic helicity 
driven case, there is a net initial total magnetic helicity injection on small scales and the dynamo 
acts to drive an inverse transfer of magnetic helicity to large scales. 

Coming back to the examples discussed in Sec. I, note that the magnetically dominated 
dynamo, or dynamical magnetic relaxation is most relevant for astrophysical coronae or laboratory 
devices, where magnetic helicity is injected into the system and the system evolves to relax 
dynamically. The kinetic helicity driven dynamo is most relevant inside astrophysical objects, 
which are not magnetically dominated. Note however that for a star or accretion disk, both 
dynamos can actually operate symbiotically: The kinetic helicity dynamo produces large scale 
fields inside the disk or star, which then rise to the corona. There the supplied fields are "small 
scale" with respect to the corona, and can subsequently dynamically relax to larger scales via 
a magnetically driven dynamo. Rapid generation of the large scale fields in this way can be an 
important source of large scale coronal magnetic fields for magnetically mediated outflows. 

Certainly our simple study herein cannot yet describe the details of astrophysical or laboratory 
applications in detail, but we have identified some unifying physical principles of magnetically 
dominated dynamos, kinetic helicity driven dynamos, and dynamical magnetic relaxation which 
can help guide further work. 

IV. Conclusions 

We have developed a nonlinear dynamical theory of magnetic relaxation, or equivalently, a 
magnetically driven large scale helical dynamo. To isolate the basic principles and to facilitate 
comparison with the simplest numerical simulations, we solved the dynamical magnetic relaxation 
equations in a closed (or periodic) system. Compared to the h v driven dynamo of Ref. [12], the 
role of h,2 and h v are reversed in dynamical magnetic relaxation: in an /12 (h v ) driven dynamo, 
subsequent growth of h v (/12) quenches the growth of h±. Also, hi grows with the same (opposite) 
sign of /12 in an /12 (h v ) driven dynamo. We have found that dynamical magnetic relaxation 
always involves an initial rapid transfer of magnetic helicity from k2 to fei, independent of Rm 
and Ry, followed by a slow, Ry dependent evolution of h\. When /12 = 1 is fixed, h\ saturates 
at h\ ~ (&2 / ki) 2 Rm / Rv ■ When /12(C)) = 1 and hi is allowed to evolve according to (18), 
hi ~ 1/2 ~ hi at maximum before resistively decaying. 

More detailed work is needed to incorporate the principles we have identified to specific 
laboratory configurations or to magnetized astrophysical coronae of stars and accretion disks. 
However, our results lead us to predict that significant magnetic relaxation will always occur on 
dynamical time scales (~ Alfven crossing time scales): The large scale magnetic helicity reaches 
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~ 1/2 the strength of the injection scale helicity during a time Tki n , after which relaxation should 
be viscously limited. The magnetic Prandtl number and the ratio of large to small scales then 
determine the ultimate saturation values. 

The theory can be fully tested with 3-D MHD numerical experiments in a periodic box in 
which magnetic helicity is injected at some relatively large wavenumber, say kf ~ 5, with an 
initially negligible velocity. The overall evolution of the magnetic helicity and kinetic helicity 
spectra can then be measured as a function of time. Future generalizations must incorporate 
boundary terms, which can be important in both astrophysical and laboratory contexts. 

EB acknowledges DOE grant DE-FG02-00ER54600. 
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Figure 1: Solutions for ki = 1, / = 1, ki = 5, (a) Rm = Rv = 200, case CI: /12 free (b) 
Rm = Rv = 2000, case CI: /12 free The curves in (a) and (b) are identified as follows: h\ 
is the thick solid line; /12 is the thin solid line; h v is the long dashed line; — Q is the short 
dashed line. 

Figure 2: (a) Same as Fig. la, but case C2: fixed hi = 1. (b) Same as Fig. lb, but case 
C2: fixed /12 = 1. (c) h\ in case C2 for late times: solid curve is for Rm = 2000, Ry = 200; 
short dashed curve is for Rm = Rv = 200; long dashed curve is for Rm = Rv = 2000. (d) 
Same as (b), but with / = 1/10. Notice that — Q now oscillates about 0. The curves in (a) 
and (b) are identified as follows: At late times (not shown) the oscillations damp, and the 
solutions are indistinguishable from the / = 1 case. 

FIGURE lab: 
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FIGURE 2ab: 



2 




r (dimensionless) 

2 , . . , 1 . . , , , . . , _ 



0} 1 ,75 - 

— I 




20 40 SO SC 100 

~L (dimensionless) 



- 13 - 




